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Numerical analysis of MHD casson fluid flow over an exponentially 
accelerated vertical plate in embedded porous medium with ramped wall 
temperature and ramped surface concentration has investigated in the current 
investigation. The flow governing dimensional velocity, temperature 
and concentration differential equations are converted into non dimensional 
form by using non-dimensional variables. Numerical solutions to 
the converted equations are obtained by finite element method. The results 
are presented graphically and in tabular form for various controlling 
parameters. In order to highlight the validity and accuracy of our present 
method, we have compared our results with the results obtained earlier. 
A very good validation of the present numerical results has been achieved. 
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I. INTRODUCTION 

Magnetohydrodynamics (MHD) is the study of the interaction between magnetic fields and electrically 
moving conducting fluids. It came into existence in the late 1930s or early 1940s. Later in nineteenth-century 
physicist Faraday conducted few isolated experiments on MHD. Afterwards astrophysicists realized how 
omnipresent magnetic fields and plasmas are throughout the universe. Then in 1942 invention Alfven waves 
created a significant attention in astrophysics. However, there was some early pioneering work by engineer 

J. Hartmann, who invented the electromagnetic pump in 1918. Hartmann also undertook a systematic theoretical 
and experimental investigation of the flow of mercury in a homogeneous magnetic field. Indeed the term 
Hartmann flow' is now used to describe duct flows in the presence of a magnetic field. Despite Hartmann's early 
researches, it was only in the early 1960s that MHD began to be exploited in engineering. However, as the interest 
in power generation declined, research into metallurgical MHD took off. Two decades later, magnetic fields are 
routinely used to heat, pump, stir and levitate liquid metals in the metallurgical industries. The key point is that 
the Lorentz force provides a non-intrusive means of controlling the flow of metals. With constant commercial 
pressure to produce cheaper, better and more consistent materials, MHD provides a unique means of exercising 
greater control over casting and refining processes. 
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The MHD flows are encountered in many geophysical, astrophysical and engineering applications. 
Hydromagnetic flows have a key role in the fields of aeronautics, stellar and planetary magnetospheres. MHD 
concepts are utilized by the engineers in the design of heat exchangers, pumps, thermal protection, in space vehicle 
propulsion control and re-entry, in creating novel power-generating systems. The purification of molten metal 
from non-metallic inclusions through the application of magnetic field is another important feature of MHD. 
All such applications of MHD give rise to investigate the problems which involves the MHD effects. 
For instance, Nadeem et al. [1] have studied MHD flow of a casson fluid over an exponentially shrinking sheet. 
Shehzad et al. [2] have analyzed the effects of mass transfer on MHD flow of casson fluid with chemical reaction 
and suction. Casson et al. [3] have investigated Effects of radiation on MHD free convection of a casson fluid 
from a horizontal circular cylinder with partial slip in non-Darcy porous medium with viscous dissipation. 
Benazir et al. [4] have discussed unsteady MHD casson fluid flow over a vertical cone and flat plate with 
non-uniform heat source/sink. Effects of fluid suction on an oscillatory MHD channel flow with heat transfer 
studied by Ahmed et al. [5]. Effects of hall current on MHD slip flow of casson nanofluid over a stretching sheet 
with zero nanoparticle mass flux studied by Aziz and Afify [6]. 

A shear thinning liquid which is assumed to have an infinite viscosity at zero rate of shear, a yield stress 
below which no flow occurs, and a zero viscosity at an infinite rate of shear is known as casson fluid. The casson 
fluid is a plastic fluid, which yields shear stress in constitutive equations. Some of the examples of casson fluid 
model are jelly, soup, honey, tomato sauce, concentrated fruit juices, drilling operations, food processing, 
metallurgy, paints, coal in water, synthetic lubricants, manufacturing of pharmaceutical products, synovial 
fluids, sewage sludge, and many others. Human blood is also considered as casson fluid. The nonlinear 
casson constitutive equation described the properties of many polymers. Furthermore, casson fluid possesses yield 
stress and has great importance in polymer processing industries and biomechanics. Hari and Harshad [7] have 
described Soret and heat generation effects on MHD casson fluid flow past an oscillating vertical plate embedded 
through porous medium. Mohyuddin and Khan [8] have studied nonlinear radiation effects on squeezing flow 
of a casson fluid between parallel disks. Akbar [9] has analyzed the influence of magnetic field on peristaltic flow 
ofa casson fluid in an asymmetric channel: application in crude oil refinement. Vajravelu et al. [10] have presented 
peristaltic transport of a casson fluid in contact with a Newtonian fluid in a circular tube with permeable 
wall. Nagarani et al. [11] have given exact analysis of unsteady convective diffusion in casson fluid flow in 
an annulus-Application to catheterized artery. A mathematical modeling for pulsatile flow of casson fluid along 
with magnetic nanoparticles in a stenosed artery under external magnetic field and body acceleration discussed by 
Priyadharshini and Ponalagusamy [12]. 

Thermal radiation is electromagnetic radiation emitted by a body by virtue of its temperature and at 
the expense of its internal energy. Thus thermal radiation is of the same nature as visible light, x-rays, 
and radio waves, the difference between them being in their wavelengths and the source of generation. 
Thermal radiation is important has applications that include high-temperature manufacturing processes 
and materials processing, improved efficiency, and more accurate design methods necessary for energy 
conversion devices, the use of new materials, hypersonic flow analysis, and others. Sajid and Hayat [13] have 
investigated the influence of thermal radiation on the boundary layer flow due to an exponentially stretching 
sheet. Bhattacharyya and Layek [14] have studied the effects of suction/blowing on steady boundary layer 
stagnation-point flow and heat transfer towards a shrinking sheet with thermal radiation. Pramanik [15] has 
explained casson fluid flow and heat transfer past an exponentially porous stretching surface in presence 
of thermal radiation. Effects of nonlinear thermal radiation on MHD chemically reacting Maxwell fluid flow 
past a linearly stretching sheet investigated by Mohan Rami Reddy et al. [16]. Yao-mathematical and 
numerical analysis of radiative heat transfer in semi-transparent media given by Chuang Han et al. [17]. 

Rate of change in temperature over time is called ramped temperature. In few investigations 
researchers considered ramped temperature profiles, the interval for ramped with respect to time in the plate 
temperature is assumed fixed i.e. 0<t'<to ( t' and to are time and characteristic time respectively) which reduces 
to 0<t<l (t being the dimensionless time) in non-dimensional form. It is to be noted that interval for ramped 
profile varies from material to material depending upon the specific heat capacity of the material. Variable 
ramped temperature profiles appear in real world situation in building air conditioning systems, fabrication 
of thin-film photovoltaic devices, phase transition in processing of materials, turbine blade heat transfer, heat 
exchangers etc. Siva Reddy Sheri et al. [18] have analyzed heat and mass transfer effects on MHD natural 
convection flow past an impulsively moving vertical plate with ramped temperature. Siva Reddy Sheri 
and Anjan Kumar [19] have discussed finite element analysis of heat and mass transfer past an impulsively 
moving vertical plate with ramped temperature. Siva Reddy Sheri et al. [20] have investigated heat and mass 
transfer effects on MHD natural convection flow past an infinite inclined plate with ramped temperature. 
Nandkeolyar et al. [21] have given exact solutions of unsteady MHD free convection in a heat absorbing 
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fluid flow past a flat plate with ramped wall temperature. Seth et al. [22] have clarified hydromagnetic 
natural convection flow with heat and mass transfer of a chemically reacting and heat absorbing fluid past 
an accelerated moving vertical plate with ramped temperature and ramped surface concentration through 
a porous medium. MHD Natural convections flow past a moving vertical plate with ramped temperature 
studied by Vijaya Bhaskar et al. [23]. 

As per authors knowledge no work has taken place to investigate the numerical analysis on MHD 
casson fluid flow over an exponentially accelerated vertical plate in embedded porous medium with ramped 
wall temperature and ramped surface concentration. In the current numerical investigation the influence 
of casson parameter and thermal radiation on hydromagnetic, incompressible, viscous electrically 
conducting, heat absorbing fluid over an exponentially accelerated vertical plate in embedded porous medium 
with ramped temperature and ramped surface concentration has been studied. The flow governing partial 
differential equations along with initial and boundary conditions are transformed into non-dimensional 
form then solved with the aid of finite element method (FEM). The effects of various pertinent parameters on 
flow pattern are displayed graphically and interpreted in detail. Verification of computations with previous 
studies in the literature is also included. The results with respect to dimensionless velocity, temperature, 
and concentration can be used as a reference for any future study. 


2. FORMULATION OF THE PROBLEM 

In Figure 1 the flow being confined to />0, where / the coordinate is is measured in the normal 
direction to the plate and v-is taken along the wall in the upward direction. Initially, at time t’= 0, both 
the fluid and the plate are uniform temperature j ' and the concentration near the plate is assumed to be cj 

at all the points respectively. At time t’>0, the plate is given an exponentially accelerated in the vertical 

direction against gravitational field with velocity an d constant heat flux • r > +T when t'<to 

^ 00 1 W 00 J . 

*0 

and TV when t’>to respectively which is there after maintained constant TV and the level of mass transfer at 
the surface of the plate is raised or lowered to , / , when t’<to and CV when t’>to respectively 

^00 + 1 C W + C 00 J , 

r o 

which is there after maintained constant CV A uniformly magnetic field of strength Bo is applied in 
the /-axis direction. Induced magnetic field is negligible. We assume that, rigid plate, incompressible flow, 
one dimensional flow, non-Newtonian fluid, free convection and viscous dissipation term in the energy 
equation is neglected. 



Figure 1. Sketch of the physical flow problem 


The constitutive equation for the casson fluid can be written as: 
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where n = e-- e-- and e-- is the (/, j)th component of the deformation rate, n is the product of the component 

of deformation rate with itself, n c is a critical value of this product based on the non-Newtonian model, fd B is 
plastic dynamic viscosity of the non-Newtonian fluid and P y is yield stress of fluid. Under above assumptions 
and taking into account the Boussinesq approximation, governing equations are given in (2-4). 


du 


' X / f 

\+-\^--ob 0 2 u’- tfy+ Pg p T (r-ri)+ Pg p c (c 1 -cy 

7 ) dy' z k ~ 


V 1 


( 2 ) 


®:.- 5 _®^: + _L£sv + «o.( r - Ti ) 

pCp dy’ 2 P°p dy' P°P 


(3) 


8C' _ d 2 C' 


hic'-c’J 


(4) 


With following initial and boundary conditions: 

u' =0,7’' = 7^, C' = C'^ for y'>0 and t’<0,u' = U 0 e a '‘, at y'=0 for f>0, 
T i = \ T ™+( T w- T ™) t ' /t o at y'=0 for 0<t'<t 0 
[K aty'=0 for t'>t 0 , 

C = 


C nn + ( C w +C co Y /t 0 at for 0 <t'<t 0 . 


C w aty'= 0 for t'>t 0 , 

u' —^ 0, T' —^ Tf , C'^Cfas y'-^oo for t'> 0 


(5) 


Introducing the following dimensionless quantities: 
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Making use of (6) in (2), (3), and (4), we get non-dimensional form as: 
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(9) 


The non-dimensional form of initial and boundary conditions are: 

u=6 = C = 0, for y > 0 and t = 0 ,u = e al 


_ ft. 

0 < t < 1) 
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3. METHOD OF SOLUTION 

The finite element method is a powerful technique for solving differential or partial differential 
equations as well as integral equations. An excellent description of finite element formulations is available in 
Bathe [24] and Reddy [25]. The basic concept is that the whole domain is divided into smaller elements 
of finite dimensions called “Finite Elements”. It is the most versatile numerical technique in modern 
engineering analysis and has been employed to study various problems in heat transfer, fluid mechanics, 
chemical processing, rigid body dynamics, solid mechanics, electrical systems and many other fields. 
The finite element method (FEM) consists of following steps: 

Step 1: Discretization of the domain 

The fundamental concept of the FEM is to divide the region of the problem into small connected 
pieces, called finite elements. The group of elements is called the finite element mesh. These finite elements 
are associated in a non-overlapping manner, such that they completely cover the entire space of the problem. 
Step 2: Invention of the element equation 

A representative element is selected from the mesh and the variational formulation of the given 
problem is created over the typical element. Over an element, an approximate solution of the variational 
problem is invented, and by surrogating this in the system, the element equations are generated. The element 
matrix, which is also known as stiffness matrix, is erected by using the element interpolation functions. 

Step 3: Assembly of the element equations 

The algebraic equations so achieved are assembled by imposing the inter element continuity 
conditions. This yields a large number of mathematical equations known as the global finite element model, 
which governs the whole domain. 

Step 4: Imposition of the boundary conditions 

On the accumulated equations, the initial and boundary conditions (10) are imposed. 

Step 5: Solution of assembled equations 

The assembled equations so obtained can be solved by any of the numerical methods, namely, Gauss 
elimination technique, LU decomposition technique, and the final matrix equation can be solved by iterative 
technique. For computational purposes, the coordinate y is varied from 0 to )W=1.5, where y max represents infinity 
i.e., external to the momentum, energy and concentration edge layers. 


4. VALIDATIONS OF NUMERICAL RESULTS 

In order to validate and verify the present numerical method used for this study in the special case is 
compared with available results of Raj Nandkeolyar et al. [21] and Seth et al. [22] is exhibited in 
Tables 1-2. These tables display the comparison of Nusselt number and Sherwood number for both ramped 
and isothermal plates and it is conveyed there is good agreement between the results exist. 


Ta ble 1. Comparison of Nusselt number (when ^ 00 , R ~ 0 and Kr - 0) 


t 

Pr 

Results by Nandkeolyar et al. [21] 

Present Results 

Ramped 

Isothermal 

Ramped 

Isothermal 

0.3 

0.71 

0.57134752 

1.11605411 

0.571347519 

1.11605410 

0.7 

0.71 

0.96929143 

0.92531051 

0.969291430 

0.92531050 

0.7 

7.00 

3.04350641 

2.90540943 

3.043506410 

2.905409429 


Table 2. Comparison of Sherwood number (when ^ 00 

and R — Q) 

t 


Results by Seth et al. [22] 

Present Results 

Kr 

Ramped 

Isothermal 

Ramped 

Isothermal 

0.3 

0.2 

0.295649 

0.525702 

0.2956500 

0.5257020 

0.5 

0.2 

0.386593 

0.428415 

0.3865929 

0.4284149 

0.3 

2.0 

0.344659 

0.839945 

0.3446600 

0.8399449 


5. RESULTS AND DISCUSSION 

Extensive numerical solutions have been found for flow governing velocity, temperature 
and concentration (6-9) under the initial and boundary conditions (10), by using finite element method. 
Obtained computational results are presented graphically in Figures 2-16 to illustrate the influence 
of governing parameters viz. Cassion parameter, Magnetic parameter, porosity parameter, thermal Grashof 
number, mass Grashof number, Prandtl number, heat absorption, thermal radiation, Schmidt number 
and chemical reaction parameter on velocity, temperature and concentration. 
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The behavior of fluid velocity against Cassion parameter y for both ramped and isothermal plates is 
illustrated in Figure 2. It reveals that an increase in y leads to diminish in velocity profile, since casson 
parameter makes the velocity boundary layer thickness shorter. It is further observed from this Figure 2 that 
plastic dynamic viscosity enlarged on increasing y it creates resistance in fluid flow and hence fluid velocity 
is decreased. Figure 3 depicts the influence of magnetic parameter M on velocity profile for both ramped 
and isothermal plates. It is perceived from this Figure 3 that velocity gets reduce by increasing M. This is due 
to the fact that the applied magnetic field normal to the flow direction induces the drag like force called 
Lorentz force which has a tendency to decelerate fluid flow. 




Figure 2. Velocity profile for casson parameter 


Figure 3. Velocity profile for magnetic 
parameter 


In Figure 4 profiles of velocity have been plotted against various values of permeability parameter K 
for both ramped and isothermal plates by keeping other parameters fixed. It is identified from this Figure 4 
that for large values of K , velocity boundary layer thickness enlarges which explains the physical situation 
that as K increases, the resistance of the porous medium is lowered which boosts up the momentum 
development of the flow regime, ultimately enhances the velocity field. The control of thermal Grashof 
number Gr on fluid velocity for both ramped and isothermal plates is demonstrated in Figure 5. Thermal 
Grashof number is the ratio of thermal buoyancy force to viscous hydro dynamic force. It is observed that 
velocity increases with increasing values of Gr. The flow is accelerated due to the enhancement in 
the buoyancy forces corresponding to the increasing values of Grashof number. 




Figure 4. Velocity profile for porosity parameter 


Figure 5. Velocity profile for Grashof number 


Int J Appl Power Eng, Vol. 9, No. 2, August 2020: 89 - 99 





















Int J Appl Power Eng 


ISSN: 2252-8792 


95 


n 


The influence of mass Grashof number on velocity profile for both ramped temperature and 
isothermal plate is displayed in Figure 6. The mass Grashof number Gm characterizes the ratio of the mass 
buoyancy force and viscous hydrodynamic force. As expected, it is observed that, for both ramped 
temperature and isothermal plate fluid velocity increases on increasing Gm. This implies that fluid velocity is 
getting accelerated due to enrichment in mass buoyancy force. Figures 7-8 portray the effect f Prandtl 
number Pr on velocity and temperature profiles for both ramped temperature and isothermal plates. 
The Prandtl Number is a dimensionless number approximating the ratio of momentum diffusivity (kinematic 
viscosity) to thermal diffusivity. 



Figure 6. Velocity profile for mass Grashof number 




Figure 7. Velocity profile for Prandtl number 


Figure 8. Temperature profile for Prandtl 
number 


From Figure 7 it is perceived that velocity profile diminishes on rising Pr. Increases in 
Prandtl number significantly depress velocity boundary layer. Since Prandtl number causes an increase in 
the viscosity of the fluid which makes the fluid thick and it leads to a decrease in the velocity. Figure 8 
reveals that temperature profile decreases on increasing Prandtl number. Physically, if Pr increases, 
the thermal diffusivity decreases and this phenomenon lead to decrease energy ability, which reduces 
the thermal boundary layer. The influence of the heat absorption Q on the velocity and temperature fields 
for both ramped temperature and isothermal plates is presented in Figures 9 and 10. The presence 
of heat absorption boundary layer generates energy which cases the temperature of the fluid to increase. 
This increasing temperature produces an increase in the flow field due to buoyancy effect. 
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Figure 9. Velocity profile for heat absorption 


Figure 10. Temperature profile for heat 
absorption 


Velocity and temperature profiles against thermal radiation for both ramped and isothermal plates 
are displayed in Figures 11-12. From Figure 11 it is seen that an increase in radiation parameter R leads to 
an increase in the fluid velocity for both ramped temperature and isothermal plates. Because when 
the strength of heat generated through thermal radiation is increased, the bond holding the components 
of the fluid particles is easily broken and the fluid velocity is increased. It is recognized from Figure 12 that 
thermal boundary layer is enhanced on increasing thermal radiation parameter R , it is due to the reason that 
thermal radiation has a tendency to enlarge the temperature profile for both ramped and isothermal plates. 




Figure 11. Velocity profile for thermal radiation 


Figure 12. Temperature profile for thermal 
radiation 


Figures 13-14 depict the influence of Schmidt number on the field of velocity and concentration. 
Schmidt number is a dimensionless number defined as the ratio of momentum diffusivity (viscosity) 
and mass diffusivity, and is used to characterize fluid flows in which there are simultaneous momentum 
and mass diffusion convection processes. As the Schmidt number increases, the concentration decreases. 
This causes the concentration buoyancy effects to decrease yielding a reduction in the fluid velocity. 
Reductions in the velocity and concentration distributions are accompanied by simultaneous reductions in 
the velocity and concentration boundary layers. 
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Figure 13. Velocity profile for Schmidt number 


Figure 14. Concentration profile for Schmidt 
number 


Effect of chemical reaction parameter on fluid velocity and concentration for both ramped 
and isothermal plates is presented in Figures 15-16. It is evident from Figure 15 that velocity profiles reduce 
on rising Kr. This is due to the fact that the chemical reaction parameter has a retarding influence on 
the solute concentration which makes a decrease in the fluid velocity due to a reduced mass buoyancy force. 
Figure 16 reveals that an increase in chemical reaction parameter causes a decrease in concentration profile. 
It may be attributed that chemical reaction parameter retards concentration boundary-layer thickness. 




Figure 15. Velocity profile for chemical reaction 


Figure 16. Concentration profile for chemical 
reaction 


6. CONCLUSIONS 

A computational solution has been carried out for investigating the thermal radiation effect on MHD 
casson fluid flow over an exponentially accelerated vertical plate in embedded porous medium with ramped 
wall temperature and ramped surface concentration. Using transformations a set of partial differential 
equations have been derived for the momentum, energy and species diffusion. Finite element method is 
engaged to solve the obtained equations. Our computational results have shown that: 

a. Velocity profile accelerated with the increasing values of porosity parameter, thermal Grashof number, 
mass Grashof number, heat absorption, and radiation parameter. On the other hand cassion parameter, 
magnetic parameter, Prandtl number, Schmidt number and chemical reaction parameter have been shown 
to retard fluid velocity for both ramped and isothermal plates. 
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b. Increasing the Prandtl number causes a deceleration in the temperature of the flow. While a reverse 
tendency has been identified against increasing values of heat absorption and radiation parameter for both 
ramped and isothermal plates. 

c. Increasing Schmidt number and chemical reaction parameter results a reduction in concentration values 
for both ramped and isothermal plates. 
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